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The Orion vehicle is being designed to provide nominal crew transport to the lunar 
transportation stack in low Earth orbit, crew abort prior during transit to the moon, and 
crew return to Earth once lunar orbit is achieved. Design of guidance and navigation 
algorithms to perform maneuvers in support of these functions is dependent on the support 
provided by navigation infrastructure, the performance of the onboard GN&C system, and 
the choice of trajectory maneuver methodology for outbound and return mission phases. 
This paper documents the preliminary integrated analyses performed by members of the 
Orion Orbit GN&C System team investigating the navigation update accuracy of a modern 
equivalent to the Apollo era ground tracking network and the expected onboard dispersion 
and navigation errors during a lunar mission using a linear covariance error analysis 
technique. 


I. Introduction 

T he Orion vehicle is being designed to provide nominal crew transport to the lunar transportation stack in low 
Earth orbit, crew abort prior during transit to the moon, and crew return to Earth once lunar orbit is achieved. 
Design of guidance and navigation algorithms to perform maneuvers in support of these functions is dependent on 
the support provided by navigation infrastructure, the performance of the onboard GN&C system, and the choice of 
trajectory maneuver methodology for outbound and return mission phases. This paper documents the preliminary 
integrated analyses performed by members of the Orion Orbit GN&C System team investigating the navigation 
update accuracy of a modern equivalent to the Apollo era ground tracking network and the expected onboard 
dispersion and navigation errors during a lunar mission using a linear covariance error analysis technique. The 
STDN navigation analysis is performed only on the transEarth trajectory to establish that a working simulation for 
Orion analysis can reproduced Apollo equivalent accuracies. This tool will then be used in support of trades for 
sizing the infrastructure needed to support the upcoming lunar exploration missions. In parallel with this navigation 
analysis, a the linear covariance study assumes Apollo ground navigation accuracies to analyze both outbound 
translunar and tranEarth return trajectories. The translunar trajectories are nominally the responsibility of 
Constellation elements other than Orion, but the dispersions from this mission phase is needed to initialize abort 
scenario development for the Orion vehicle. 

The CEV/LLV stack will be inserted into a translunar orbit by the Earth Departure Stage. Soon after translunar 
injection (TLI) maneuver the EDS will separate from the CEV/LLV stack. It will perform 4 TCMs during the 3 -day 
coast. Depending on the desired lunar orbit, a sequence of (up to three) maneuvers will be performed to insert the 
CEV/LLV stack into a low lunar orbit, that passes over the landing site. This insertion maneuver sequence could 
take up to 24 hours. After insertion, the LLV will separate from the now un-crewed CEV, which will continue 
unmanned operations until the LLV ascent stage returns to dock with the CEV. After the crew and cargo is 
transferred to the CEV, the LLV will separate from the CEV and the CEV will prepare for the Trans-Earth Injection 
(TEI) sequence. This could be a three-burn sequence with the transfer orbit between the first and the last maneuver 
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allowing for a period of upto 48 hours. After the TEI sequence, the remaining transfer time would be 72 hours to 
reach the Entry Interface point. 

The paper is organized as follows. Section 2 contains a description of the tools used in the design and analysis as 
well as the assumptions made. Section 3 contains the translunar design and analysis, including the TLI and lunar 
orbit insertion phase. Section 4 contains the design and analysis of the transEarth phase, including the TEI 
maneuvers. Finally, Section 5 contains some concluding comments. 

II. Tools and Assumptions 
A. Linear Covariance Analysis and Tools 

Linear covariance analysis is a methodology to obtain the navigation performance (statistics) of a system in a single 
run! It is a complement to Monte Carlo analysis in that it provides the statistical behavior of a GNC system without 
the thousands of Monte Carlo runs needed in order to truly represent the statistics of a complex dynamical system. 
During the early phase of a system design it is invaluable in being able to provide the engineer with the system 
performance, especially as it relates to navigation and guidance issues. It is used to quickly (and accurately) 
understand how various navigation sensors affect the navigation accuracy of a particular system design under 
consideration. Additionally, it also provides trajectory dispersion information (the difference between the actual 
trajectory and the nominal trajectory) along with variation in fuel (or AT) usage. Finally, it is the perfect tool in 
order to trade the performance of different sensors or to determine the contribution of a particular sensor (error) 
parameter. 

However, it is necessary to keep in mind the limitations of linear covariance analysis, particularly as it relates to 
rendezvous analysis. As the name suggests, it is linear covariance analysis. It relies on the methodology of Kalman 
filtering (as was originally obtained by Kalman in 1960); it operates on perturbation states, which are assumed to be 
linear. The nominal dynamics are propagated using the nonlinear dynamics. The measurement sensitivities are the 
partial derivatives of the nonlinear measurements with respect to the states of interest (as in a Kalman filter). As the 
second part of the name denotes, this analysis is based upon the covariance matrices associated with the Kalman 
filter. Recall that within virtually all variants of Kalman filters, the covariance matrices are propagated in a linear 
fashion (either by means of a state transition matrix or by integrating the Ricatti differential equation). We choose 
the latter formulation for propagating the covariance matrix. The nominal states are propagated using the full 
nonlinear equations. 

In the cislunar linear covariance analysis we deal with the trajectory dispersion covariance (the covariance of the 
difference between the actual states and the nominal states), the navigation dispersion covariance (covariance of the 
difference between the navigated states and the nominal states) and the on-board navigation errors (covariance of the 
difference between the navigation states and actual states). We carry along and operate on all these three covariance 
matrices. 

7. Dynamics Models 

For this CEV cislunar analsysis, the gravity fields of the bodies (Sun, Earth and Moon) were assumed to be 
spherical. However, all three bodies were assumed to be active during the entire trajectory. The integrator used for 
the propagation of the covariance dynamics was a fourth-order Runge-Kutta method. The trajectory of the CEV was 
propagated using a (fourth-order) Encke-Nystrom integrator. In order to retain precision in the covariance 
propagation, a 10 second time step was used. 

The attitude was modeled as a right-handed quaternion (with scalar first). For the covariance associated with the 
attitude dynamics, the Bortz formulation was used. Unlike the quaternion representation in a Kalman filter, it 
provides for an attitude covariance matrix of full rank without any approximations. 

For this early analysis, the attitude control system was not modeled. Once again, this can be incorporated as the 
system matures. Particularly for this phase of the mission, the attitude control was not considered to be particularly 
important when PTC (Passive Thermal Control, a.k.a. the barbeque roll) will, most likely, be used. 

One of the significant factors in the translunar trajectory is the translation noise. This manifests itself in 
imperfectly coupled rotational (attitude) maneuvers, venting, etc. In the filter this is modeled as process noise and is 
one of the parameters that is varied. There are two types of periods modeled: active and quiescence. For the 
nominal case, the active period was modeled as 4.5 jig(Earth) (0.027 ft/sMir). The quiescent period was modeled 
with a process noise of 0.45 qg(Earth) (0.0027 ft/sMir). The active periods last for 16 hours per day and the 
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quiescent periods last for 8. The first quiescent period begins 9 hours after TLI and these alternating 
active/quiescence periods continue for the remainder of the transfer. 

There are 12 states associated with the CEV dynamics (3 position, 3 velocity, 3 attitude, 3 attitude rate). The 
maneuvers errors are modeled as biases representing thruster misalignments (0.01 degrees (la)), scale factor (10 
ppm § (la)), and biases (0.00328 ft/s (la)), along with a random noise of 0.00328 ft/s (la). 

It is assumed that the TLI maneuver will have active guidance so that the maneuver errors and cut-off errors will 
be minimal. As a point of comparison, during Apollo, the S-IVB TLI burn was so accurate that TCM 1 was not 
performed. It is expected that the residuals associated with this maneuver would be on the order of 1 m/s. 


2. Sensor Models 

The CEV will have an extensive suite of navigation instruments. While most of them are used for rendezvous, 
the IMU, GPS receiver, and star tracker are used for all phases of flight, including cislunar operations. 

The gyros were modeled after the GG1320 RLGs in the Honeywell MIMU. It has a long-term bias stability of 
0.12 deg/hr (3a) and a random walk of 0.025 deg/Vhr. There were 18 states associated with the gyros that were 
modeled. They were: gyro drift Markov (3 states), gyro drift bias (3 states), gyro scale factor Markov (3 states), 
gyro scale factor bias (3 states), gyro misalignments (6 states). The Markov states had time constants of 1 hour and 
strength corresponding to the bias states. 

The CEV GPS receiver performance was modeled after a Viceroy class GPS receiver. It was assumed to provide 
deterministic GPS measurements accurate to 50 ft and 0.067 ft/s (la) with updates every 10 seconds. In addition to 
the measurement nose, the GPS receiver was modeled in the onboard system using 6 Markov states to model the 
GPS position and velocity errors. The maximum altitude of the GPS receiver was conservatively chosen to be 6475 
n.m. (12,000 km). 

The star tracker was modeled after a Goodrich HD 1 003 and provides attitude updates in the form of an inertial- 
to-sensor quaternion every 10 minutes. The measurement accuracy was 62 arc-seconds (3a). The misalignment 
was 36 arc-seconds (3a). 

There are infrequent updates to the navigation states of the CEV. These are performed prior to each maneuver. 
The accuracy of the ground updates is listed in Table 3, reflective of ground tracking capabilities during the Apollo 
10 mission . It is expected that today’s capabilities are comparable due to the fact that there are few tracking 
stations available and the geometry of those is not as favorable; however technology has advanced in the intervening 
years. Table 1 contains the ground-tracking accuracies used in this analysis. In upcoming years, further analysis 
will be performed to better characterize the expected cislunar navigation accuracies. 


1 a Error 

TCM 

TCM 

TCM 

TCM 

LOI 

LOI 

LOI 

Component 

1 

2 

3 

4 

1 

2 

3 

Radial (U) 

482 ft 

1502 ft 

2832 ft 

4031 ft 

3032 ft 

516 ft 

3032 ft 

Downrange (V) 

7263 ft 

13549 ft 

19349 ft 

4651 ft 

5051 ft 

18002 ft 

5051 ft 

Crossrange (W) 

13919 ft 

10020 ft 

14228 ft 

18960 ft 

8109 ft 

45571 ft 

8109 ft 

Radial Rate(U-dot) 

0.058 ft/s 

0.82 ft/s 

0.015 ft/s 

0.027 ft/s 

3.241 ft/s 

0.014 ft/s 

3.241 ft/s 

Downrange Rate (V-dot) 

0.214 ft/s 

0.39 ft/s 

0.064 ft/s 

0.008 ft/s 

1.177 ft/s 

0.343 ft/s 

1.177 ft/s 

Crossrange Rate (W-dot) 

0.318 ft/s 

0.49 ft/s 

0.061 ft/s 

0.042 ft/s 

7.300 ft/s 

1.010 ft/s 

7.300 ft/s 


Table 1: Translunar Ground Update Navigation Accuracy Parameters 
B. Design Reference Mission Trajectories 

The translunar trajectory was designed to deliver the CEV to a low lunar orbit five days after the TLI (Trans- 
Lunar Injection). There is a three-burn Lunar Orbit Insertion (LOI) maneuver sequence to insert the CEV/LSAM 
stack into a near-polar Low Lunar Orbit (LLO). For the return trip, a three-burn TEI sequence essentially unwinds 
the LOI sequence and inserts the CEV into a departure hyperbola that results in a desired Entry Interface Point. The 
maneuver times, their objectives and the nominal A Vs are given in Table 2. The corresponding transEarth (Moon- 
to-Earth) maneuver information is listed in Table 3, with the times referenced to the TEf maneuver 


§ Parts per million 

Lewis, L.R., et.al., Apollo 10 Spacecraft Dispersion Analysis, Volume 1- Navigation Error Analysis, MSC 
Internal Note No. 69-FM-83, April 14, 1969. 
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Maneuver 

Type 

Maneuver Time 
(hrs) 

Maneuver Objective 

Nominal A V (ft/s) 

TLI 

0 

Target for h p . Energy, Inclination at 
Pericynthion 

10,308 

TCMj 

6 hours (TLI + 6 hours) 

Target for h p , Energy, Inclination at 
Pericynthion 

0.0 

tcm 2 

24 hours (TLI + 1 day) 

Target for h p . Energy, Inclination at 
Pericynthion 

0.0 

tcm 3 

72 hours (LOI - 1 day) 

Target for h p . Energy, Inclination at 
Pericynthion 

0.0 

tcm 4 

90 hours (LOI - 6 
hours) 

Target for h p . Energy, Inclination at 
Pericynthion 

0.0 

LOIj 

96 hours 

Apocynthion Altitude 

826 

LOI 2 

99.5 hours 

Wedge Angle, Energy, Pericynthion 
altitude 

216 

LOI 3 

120 hours 

Apocynthion Altitude 
(circularization) 

1856 


Table 2: Translunar Maneuver Information 


Maneuver 

Type 

Maneuver Time (hrs) 
(Past TEI0 

Maneuver Objective 

Nominal A V (ft/s) 

TEIj 

2.68 

Raise Apocynthion altitude 

6147 

tei 2 

15.16 

Wedge Angle Change, 
Energy Adjust, Pericynthion 
Altitude Adjust 

1529 

tei 3 

24.05 

Target for Entry FPA, 
Azimuth, Latitude, Longitude 

3623 

tcm 5 

42.05 hours (TEI + 18 hours) 

Target for Entry FPA, 
Azimuth, Latitude 

0.0 

tcm 6 

92.03 hours (El - 16 hours) 

Target for Entry FPA, 
Azimuth, Latitude 

0.0 

tcm 7 

103.03 hours (El - 5 hours) 

Target for Entry FPA, 
Azimuth, Latitude 

0.0 


Table 3: TransEarth Maneuver Information 


III. Translunar Navigation and Dispersion Analysis 

In order to analyze the overall system performance, several trades were performed to ascertain the sensitivity to 
factors such as frequent ground updates, maneuver execution errors and vehicle noise. Whereas analysis has been 
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performed for various ground update frequencies, the best performance is obtained, as expected, for the case of 
frequent ground updates. Figure 1 shows the position and velocity navigation errors for the out-bound trajectory. 
Noted on the plots is the effect of ground navigation updates on the on-board navigation errors. 




Figure 1: On-board Position and Velocity Navigation Errors 

7. Navigation Error Analysis 

The UVW coordinate system, while useful in certain applications (such as rendezvous) are less useful for orbits 
with high eccentricity (such as the Earth-Moon and Moon-Earth transfers). As such it is useful to find other metrics 
for analyzing navigation (and dispersion) errors. One of the most useful metrics is to propagate the (current) 
navigation errors to the final epoch (usually periapsis or entry- interface). For covariance analysis, this is done by 
means of the state transition matrix and a transformation from inertial coordinates to the parameters of interest. For 
the translunar trajectory, the quantities of interest are the pericynthion (closest approach) altitude, inclination, and 
energy. Figure 2 contains a plot of the navigation error mapped to pericynthion altitude, with the right plot 
containing the zoomed in mapped pericynthion altitude error. Figure 3 contains the mapped pericynthion inclination 
error and the Cj error ( C 3 = —/d/2 a). 


CEV 3-o On-Board Navigation Errors Mapped to Pericynthian Altitude 




Figure 2: Translunar coast navigation errors mapped to pericynthion altitude 
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CEV On-Board 3-o Inclination Navigation Error 


CEV 3-o On-Board Navigation Errors Mapped to C 3 at Pericynthion 




Figure 3: Translunar coast navigation errors mapped to Pericynthion Inclination and C 3 

Figure 4 shows the periapsis altitude and semi-major axis navigation error after the first LOI maneuver, when the 
CEV is captured into lunar orbit. Figure 5 shows the inclination and longitude of the ascending node navigation 
error upon capture in lunar orbit. 


CEV On-Board 3-o Pericynthian Altitude Navigation Error CEV On-Board 3-o Semi-major Axis Navigation Error 




Figure 4: Periapsis alitude and Semi-major axis navigation errors 


CEV On-Board 3-o Inclination Navigation Error CEV On-Board 3-o Longitude of the Ascending Node Navigation Error 




Figure 5: Inclination and Longitude of the Ascending Node navigation errors 
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2. Dispersion Analysis 

Linear covariance analysis also provides the analyst with trajectory dispersions (defined as the difference 
between the actual trajectory and the nominal trajectory). As with the navigation errors, the translunar trajectory 
dispersion errors were mapped to pericynthion and expressed in terms of the quantities of interest at pericynthion. 



Time (hrs) 


CEV 3-o Trajectory Dispersions Mapped to Pericynthian Altitude 



Figure 6: Translunar Trajectory Dispersions Mapped to Pericynthion Altitude Errors 


Figures 8 and 9 contain the trajectory dispersions mapped to orbital quantities of the post-insertion orbit. One can 
see that the dispersions after the sequence of insertion maneuvers are completed are quite good and well-contained. 
This bodes well for the translunar and insertion phase of the mission, assuming that the ground navigation system is 
kept to the accuracy of the Apollo tracking station network. Finally, Table 4 contains the AV dispersion statistics 
associated with the maneuvers along this trajectory and Table 5 contains the final post- insertion orbit. 


CEV 3-0 Inclination Dispersion Error x 10 < CEV 3 "° Position Dispersion Mapped to C 3 at Pericynthion 




Figure 7: Translunar Trajectory Dispersions Mapped to Pericynthion Inclination and C3 Errors 
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CEV 3-o Pericynthian Altitude Dispersion Error (nautical miles) 


2.6 


CEV 3-o Pericynthian Altitude Dispersion Error 


CEV 3-o Semi-major Axis Dispersion Error 




Figure 8: Post-insertion Trajectory Dispersions Mapped to Semi-major axis and Pericynthion Altitude 

CEV 3-o Inclination Dispersion Error CEV 3-o Longitude of the Ascending Node Dispersion Error 




Figure 9: Post-Insertion Trajectory Dispersions Mapped to Inclination and Longitude of the Ascending Node 


Maneuver Type 

Nominal AV 
(ft/s) 

°AF 

(ft/s) 

TCMi 

0 

0.997 

tcm 2 

0 

0.140 

tcm 3 

0 

0.496 

tcm 4 

0 

0.423 

LOIi 

825.997 

9.360 

loi 2 

215.725 

3.485 

loi 3 

1855.804 

61.297 


Table 4: la AV Dispersions 


Error 

Component 

(la) 

Semi-major Axis 

677.2 ft 

Periapsis Altitude 

1426.5 ft 

Inclination 

0.0229 deg 

Longitude of the Ascending Node 

0.0206 deg 


Table 5: la Dispersions at the Conclusion of the LOI phase (t LO i 3 + 1 hour) 
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IV. Trans-Earth Navigation and Dispersion Analysis 


As in the case of the translunar analysis, both the navigation and the trajectory dispersions will be analyzed. The 
navigation errors for the entire phase (in terms of UVW) are illustrated in Figure 10. 


CEV 3-o Position Navigation Error in UVW CEV 3-o Velocity Navigation Error in UVW 




Figure 10: TransEarth On-board Navigation Errors 

7. Navigation Error Analysis 

As in the previous section, the navigation errors will be analyzed as a function of the sub-phase: lunar orbit and 
trans-Earth. This time, however, we look at the TransEarth Injection phase (which includes the TransEarth Injection 
(TEI) maneuvers) followed by the TransEarth coast phase. Figure 1 1 contains the plots of the pericynthion altitude 
and semi-major axis navigation errors; Figure 12 contains the inclination and longitude of the ascending node 
navigation errors. One can clearly see the effect of the ground (Earth-based) navigation updates on the navigation 
state. Figures 13 and 14 contain the on-board navigation errors mapped to the entry- interface point (400,000 ft 
above the surface of the earth). One can see the effect of the ground-based navigation updates and the effect of GPS 
acquisition. For this analysis, it was assumed that GPS was acquired l A way through the GPS shell (-12,000 km 
above the surface). 


CEV On-Board 3-o Pericynthian Altitude Navigation Error 



CEV On-Board 3-cr Semi-major Axis Navigation Error 



0 5 10 15 20 25 30 

Time (hrs) 


Figure 11: Pericynthion Altitude and Semi-major axis navigation errors for the Trans-Earth Injection Phase 
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CEV 3-o Inclination Navigation Error (degrees) 


CEV On-Board 3-o Inclination Navigation Error 


CEV On-Board 3-o Longitude of the Ascending Node Navigation Error 



Time (hrs) Time (hrs) 


Figure 12 : Inclination and longitude of the ascending node navigation errors for the TEI Phase 




Figure 13: TranslEarth coast navigation errors mapped to entry interface flight path angle and heading angle 

CEV 3-o On-Board Navigation Errors Mapped to Entry Interface Latitude CEV 3-o On-Board Navigation Errors Mapped to Entry Interface Longitude 




Figure 14: TransEarth coast navigation errors mapped to entry interface latitude and longitude 

2. TransEarth Dispersion Analysis 

The trajectory dispersions for the TEI maneuver phase (the phase between TEIi and TEI 3 ) of the pericynthion 
altitude and semi-major axis errors are shown in Figure 15 and the inclination and longitude of the ascending node 
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dispersions are illustrated in Figure 16. For the transEarth coast phase, Figure 17 shows the trajectory dispersions 
mapped to the entry interface flight-path angle and heading angle and Figure 18 shows the trajectory dispersions 
mapped to entry latitude and longitude. 

Finally, Table 6 has the entry-interface trajectory dispersions expressed in terms of entry flight path angle, 
heading angle, latitude and longitude and Table 7 has the Is DV statistics. One can see that the dispersions at TEI3 
are the largest, but that is anticipated due to the large size of that maneuver. It is primarily caused by navigation 
errors during the prior maneuvers. The TCM dispersions for the transEarth coast phase are well contained. 


CEV 3-o Pericynthian Altitude Dispersion Error CEV 3-o Semi-major Axis Dispersion Error 




Figure 15: Pericynthion altitude and semi-major axis dispersion errors for the TEI maneuver phase 



Time (hrs) 


CEV 3-o Longitude of the Ascending Node Dispersion Error 



Figure 16: Inclination and longitude of the ascending node dispersion errors for the TEI phase 
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Figure 17: TransEarth coast trajectory dispersion errors mapped to entry interface flight path angle and 

heading angle 


CEV 3-o On-Board Navigation Errors Mapped to Entry Interface Latitude CEV 3-o On-Board Navigation Errors Mapped to Entry Interface Longitude 



Time (hrs) Time (hrs) 


Figure 18: TransEarth coast trajectory dispersion errors mapped to entry interface latitude and longitude 


Error 

Entry Dispersion 

Component 

(la) 

Flight Path Angle 

0.0486 deg 

Heading Angle 

0.0234 deg 

Latitude 

0.0758 deg 

Longitude 

0.2089 deg 


Table 6: lo Dispersions at Entry Interface 


Maneuver Type 

Nominal AV 
(ft/s) 

(ft/s) 

TEli 

1873.634 

8.376 

tei 2 

465.996 

8.647 

tei 3 

1104.216 

25.165 

tcm 5 

0 

0.713 

tcm 6 

0 

0.869 

tcm 7 

0 

0.442 


Table 7: lo AV Dispersions 
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V. Conclusion 


This paper has documented the analysis and design performed for the Orion translunar and transEarth phases. 
The analysis has demonstrated that if the Apollo-era tracking and data network navigation data is available, 
excellent AV dispersions and navigation performance is obtained. Future analysis will further flesh out the details 
as the design matures. 
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